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1. Introduction 



Deep inelastic structure functions of mesons and nucleons are interesting since those functions 
give us information about the momentum and spin carried by quarks and gluons inside hadrons. 
Lattice QCD can calculate moments of structure functions from first principles. By doing an op- 
erator product expansion on the hadronic tensor, the moments of structure functions are related to 
reduced matrix elements of certain local operators. For example, for the case of scattering of a spin 
one particle, the details are given in Ref.[^. On the lattice, calculations on moments of structure 
functions of mesons and nucleons can be found in, e.g., Ref.[Q |3[, Q ^]. 

The Wilson twisted mass formulation of lattice QCD was introduced in Ref.[§, ^. It provides 
automatic &{a) improvement when tuned to maximal twist, which can be achieved by setting the 
PCAC mass to zero. Also, twisted mass fermions are protected from unphysical fermion zero 
modes, thus the problem of exceptional configurations is avoided. Quenched simulations with Wil- 
son twisted mass fermions have been shown to be successful||8|, |, 1^, 11]. Specifically, a quenched 
calculation of the lowest moment of the quark distribution function (x) in a pion was done in 



Ref.||12|]. Recently, two flavor dynamical simulations []13|] with Wilson twisted mass fermions gave 
precise results on low energy constants in the chiral effective Lagrangian. Using those configura- 
tions, in this work we calculate (x) in a pion. The computation is performed at pion masses in the 
range of ~ 300 — 500 MeV. In Table [l|, the parameters of our simulations are collected. The gauge 
action used in the simulations is the tree-level Symanzik improved gauge action. For more details 



of our data, please see Ref.[|14|]. 





aCfm) ajJL 


m;j(GeV) 


L^xT 


A'meas 


(x)bare 


3.9 


0.0855(6) 0.0100 


0.4839(12) 


24^ X 48 


170 


0.295(3) 




0.0085 


0.4470(12) 


24^ X 48 


167 


0.279(4) 




0.0064 


0.3903(9) 


24^ X 48 


250 


0.279(5) 




0.0040 


0.3131(16) 


24^ X 48 


230 


0.268(8) 




0.0040 


0.3082(55) 


323 X 64 


316 


0.251(8) 



Table 1: Simulation parameters, the number of measurements(A'meas) and preliminary results of (x)^^^^. 
Measurements were done on gauge configurations separated by 20 HMC ttajectories with trajectory length 



2. Methodology 

By using the optical theorem, the cross section of deep inelastic scatterings can be related to 
the hadronic tensor or the imaginary part of the forward current-hadron scattering amplitude W^^, 
which is given by 

W^'''{p,q,X,X') = ^J d'xe'^-^p,X'\[f{x),fm\p,X). (2.1) 

Here A and A' are the polarization of the target hadron. p and q are the momenta of the hadron and 
the virtual photon. Depending on the spin of the hadron, the above tensor W^^ can be decomposed 
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into a number of independent structure functions. The spin-averaged structure functions F\{x,Q^) 
and F2{x,Q^) can tell us overall densities of quarks and gluons in a hadron. In the parton model, 
to leading order, the single flavor structure function F\ (x) is half the probability of finding a quark 
with momentum fraction x. If defining the «th moment of a function f{x) as 

Mn{f) = f'x"-'f{x)dx, (2.2) 
Jo 

then for the pion, to leading twist order, we have from the operator product expansion of W^^ 

2Mn{Fi)=d^Kn, Mn-i{F2)=dn\n, (?i>2,even), (2.3) 

where C^n^ = 1 + &{ois) are the Wilson coefficients, and the reduced matrix element v„ is defined 
by 

(^l^^^'"''"^ -traces|p) = 2vn[p^^ ■ ■ ■ /j''" - traces] . (2.4) 
Here { • • ■ } means symmetrization on the Lorentz indices and the twist-2 operators are given by 

^n,-ii„ = (^)n-i^^.^,^^,^/.. . . (2.5) 

where D=D — D and Gff is a diagonal flavor matrix. 

To compute the lowest moment of the quark distribution function (a), we need to consider 
« = 2 in Eq.(2.3). The corresponding twist-2 operator has two irreducible representations on the 
lattice[15]. In our calculation we use the following operator 

1_ 1 

^44(x) = -u{x)[y4 Da -t £ Yk DkHx), (2.6) 
^ ^ k=i 

where = ^(Vm + V/j) with (Vp being the usual forward (backward) derivative on the 
lattice. With the above operator, no external momentum is needed in our calculation, which is 
advantageous since an external momentum increases the noise to signal ratio. One can also use 
an operator from the other representation, which needs an external momentum. In the continuum 
limit, the two operators should give the same result on (x) . 
From Eq.(2.4) the bare moment (x)^^^ is given by 

(^)bare ^ ^ _L^^^o\ff44\7i,0), (2.7) 

where the matrix element (7r,0|(^44|;r,()) between two pions at rest is calculated from the ratio of 
the following 3-point function and 2-point function with a source at f = and a sink atT/2: 

{n,d\ff44\n,d)=4m^-^^^ (0«f«r/2). (2.8) 

Here 

C44it)='£{PSiT/2,x)^u{t,y)PS\0,0)), (2.9) 

xy 

C^{T/2)='£{PSiT/2,x)PS\0,0)), (2.10) 
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and PS{x) = u{x)y^d{x) is the interpolating field for the pion. There are two contributions in the 
Wick contractions of C44(f): one connected diagram and one disconnected diagram. The discon- 
nected contribution is ignored in our calculation at this moment, but will be computed by using a 
stochastic source method. 

The above two and three point correlators are evaluated by using a stochastic time slice source 



(Z(2)-noise in both real and imaginary part) [ ]16| , |17| , |18| ] for all color, spin and spatial indices, i.e., 
the quark propagator X^{y) is obtained by solving 

Y^D^^^^{z,y)xl{y) = ^{z)l5,^_o (source at ? = 0), (2.11) 

where the Z(2) random source E, {z)a satisfies the random average condition 

{^*{xU{yfp) = 5,s5aM5a^p. (2.12) 
The generalized propagator| 19] ^^(y) needed in computing €44(1) is obtained by solving 

X<^(z,3')rJ(3') = 75^^(z)5,o,r/2 (sink at t = T/2). (2.13) 
y 

The advantage of using the above stochastic source compared with a point source is that much 
less inversions are needed. With a point source, 24 inversions per gauge configuration are needed: 
12 (3 colors X 4 spins) for the quark propagator and 12 for the generalized propagator With the 
stochastic source, only two inversions are needed. One for the quark propagator, another one for 
the generalized propagator. Note that this saving ratio in the number of inversions is special for the 
pseudoscalar meson. For other meson correlators, smaller saving ratios can be achieved. 

In our data analysis, measurements are done on gauge configurations separated by 20 HMC tra- 



jectories with trajectory length T = 1/2. Statistical errors are from a Gamma-function analysis[2C]. 



3. Preliminary results and outlook 

At the smallest quark mass, ajj. = 0.0040, we calculated (x) using both a point source method 
and the above stochastic source method on the 24^ x 48 lattice. The results are in agreement within 
errors. Figjl] shows the comparison of the pion effective masses obtained from the two methods 
on 230 gauge configurations. As we can see, the statistical errors from the two methods are of the 
same size. In Fig. we compare the results of (x)^^^ obtained by using the two sources. The two 
plateaus of (x)^^^ around T /4 and 3r/4 are averaged to increase statistics. Again, we find the 
statistical errors from the two methods are comparable. However, the stochastic source method is 
much cheaper in computer time. 

Our preliminary results of (x)^^^ are collected in the last column of Table |l|. We only give 
bare quantities here since the renormalization constant of the matrix element has not been calcu- 
lated yet, which we plan to compute non-perturbatively. Fig. |3| shows the results of (x)^^^ against 
the pion mass squared. The cross in the graph at the lowest quark mass is from the 32^ x 64 lat- 
tice. It shows that the finite lattice volume effects are not big. A linear extrapolation in gives 
bare _Q 246(10) in the chiral limit (using only the points from the 24^ x 48 lattice). 
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Figure 1: Comparison of the pion effective masses for a/i = 0.004 obtained from a point source and a 
stochastic source on the 24^ x 48 lattice. They are from the same 230 gauge configurations. 
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Figure 2: {x)^^^ for api — 0.004 obtained from a point source and a stochastic source. We show the average 
of the two plateaux around T /A and 3r/4. The stochastic source method is much cheaper in computer time. 



Using Wilson twisted mass fermions, our two flavor dynamical simulations go down to a 
pion mass of around 300 MeV. Our first attempt in this work shows we can get (a) with small 
statistical errors of a few percent. This will enable us to compare numerical results of (x) in 



pions and nucleons with predictions from the chiral perturbation theory (see for example ||21|, |22|, 
23]). With simulations at two more lattice spacings, we will be able to do an extrapolation to 
the continuum limit (the data analysis is in progress). For the future, we plan to compute the 
disconnected diagrams using a stochastic source method and calculate the renormalization constant 
non-perturbatively. We also plan to compute (x) in nucleons. 
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Figure 3: {x) against the pion mass squared. Finite size effects at the lowest pion mass are not big. 
A linear extrapolation gives (.ji:)^^^=0. 246(10) in the chiral limit (using only the points from the 24^ x 48 
lattice). 
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